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Abstract 

This  paper  considers  subspace  based  adaptive  detection 
in  the  context  of  the  likelihood  ratio  test  studied  by  Kelly, 
[  1 }.  The  probability  of  false  alarm  for  this  test  depends  only 
on  the  subspace  dimension  while  the  probability  of 
detection  is  a  function  of  the  subspace.  Thus,  we  propose 
choosing  the  transformation  onto  the  subspace  to  maximize 
the  probability  of  detection  over  a  likely  class  of  noise  and 
interference  scenarios.  An  approximate  solution  to  this 
optimization  problem  is  described.  This  approach  can  lead 
to  dramatic  increases  in  the  probability  of  detection  given  a 
fixed  number  of  riaia  observations  due  to  a  large  gain  in  the 
statistical  stability  associated  with  the  reduced  dimension 
subspace.  The  relationship  between  subspace  design  for 
adaptive  detection  and  partially  adaptive  beam  former  design 
is  explored.  Simulations  verify  the  analysis. 

1.  Introduction 

Processing  of  sensor  array  data  in  subspaces, 
Sometimes  termed  the  beamspace  domain,  has  been 
proposed  for  adaptive  beamforming  and  source  location 
estimation  (e.g.  [2-4]).  Benefits  of  mapping  the  data  into  a 
subspace  prior  to  processing  include  reduced  real  lime 
computational  burden,  improved  adaptive  convergence 
behavior,  increased  resolution  and  robustness.  In  this  paper 
we  consider  subspace  based  adaptive  detection  for  sensor 
array  data  in  the  context  of  the  likelihood  ratio  test  studied 
by  Kelly  [1].  The  signal  to  be  detected  is  modeled  at  the 
sensor  outputs  as  a  known  vector  with  unknown  amplitude. 
The  noise  and  interference  are  modeled  as  zero  mean 
Gaussian  random  vectors  with  unknown  covariance  matrix. 
The  likelihood  ratio  test  is  dependent  on  the  unknown  signal 
amplitude  and  noise/interference  covariance,  so  the  detection 
statistic  is  obtained  by  maximizing  the  likelihood  ratio  over 
these  unknown  parameters  given  the  data.  Kelly  [1]  shows 
that  the  false  alarm  rate  for  this  statistic  is  independent  of 
the  noise/interference  covariance  matrix  and  derives  an 
expression  for  the  probability  of  detection. 

Kelly's  [1]  results  are  generalized  here  for  detection  in 
arbitrary  reduced  dimension  subspaces.  Expressions  for  the 
probabilities  of  detection  and  false  alarm  as  a  function  of  the 
subspace  are  given.  The  probability  of  false  alarm  is 
dependent  only  on  the  subspace  dimension  so  we  propose 
choosing  the  subspace  to  maximize  the  probability  of 
detection  over  a  likely  class  of  interference  scenarios.  An 
approximate  solution  to  this  optimization  problem  is 
obtained  by  choosing  the  components  of  the  subspace  one  at 
a  time,  with  each  additional  component  chosen  to  maximize 
the  probability  of  detection.  In  spite  of  the  loss  in  SNR 
which  results  from  subspace  processing,  dramatic 
improvements  in  detection  performance  can  be  obtained 
because  of  the  increased  statistical  stability  associated  with 


estimation  of  the  covariance  matrix  in  the  reduced 
dimensioned  subspace. 

The  outline  of  this  paper  is  as  follows.  Section  2 
introduces  the  detection  problem  and  modifies  Kelly's  results 
for  subspace  processing.  The  relationship  between  partially 
adaptive  beamforming  and  the  parameters  of  the  subspace 
detection  problem  is  also  illustrated.  A  technique  for 
choosing  the  subspace  to  approximately  maximize  the 
probability  of  detection  over  a  set  of  likely  interference 
scenarios  is  described  and  analyzed  in  section  3.  Several 
simulation  examples  are  given  in  section  4  to  illustrate  the 
effectiveness  of  this  approach.  A  summary  is  provided  in 
Section  5. 


2.  Adaptive  Detection 

Let  the  N  dimensional  vector  xj  be  the  data  vector 
observed  at  the  N  sensor  outputs  at  lime  i.  The  x;  are 
assumed  to  be  independent  and  Gaussian  distributed  with 
means  under  the  two  hypotheses 

Ho:  E{x.}  =  0,  i  =  0.1 . K 

H,:  Ejxoj-bs.EfxJ-O,  U . K 

The  covariance  is  the  same  under  both  hypotheses 
H0.H,:  e{(x.  -  Ejx^Jfx,  -  E{x. ))”}  *  M 

The  signal  vector  s  is  assumed  known,  but  the  signal 
amplitude  parameter  b  and  the  noise  covariance  M  are 
unknown.  Given  K  signal  free  data  vectors  xj,  i  =  1.2....K, 
we  wish  to  determine  whether  or  not  a  signal  is  present  in 
the  data  vector  xo- 

Kelly  [1]  uses  a  likelihood  ratio  approach  wherein  the 
unknown  parameters  M  and  b  are  replaced  by  their 
maximum  likelihood  estimates.  This  leads  to  a  threshold 
based  decision  statistic  whose  probability  of  false  alarm 
(PFA)  is  dependent  only  on  K,  N  and  the  threshold.  An 
expression  for  the  probability  of  detection  is  also  derived  in 
[lj.  This  approach  is  repeated  below  assuming  the  data  is 
transformed  into  a  subspace  prior  to  performing  detection. 

Define  an  N  by  P  (P<  N)  dimensioned  full  rank  matrix 
T  and  transformed  data  vectors  zj  =  THXj,  i  *  0,1 . K. 


(1) 

(2) 


The  subspace  data  remains  zero  mean  except  for  zo  under  Hj 
which  has  mean  E{zo)  *  bTHs  *  bsj.  The  transformed 
covariance  matrix  is  My  =  T^MT.  The  joint  probability 
density  functions  fp  and  f]  under  Hq  and  Hi  are  given  by 


sNdetM,. 


«xp[-c{M^Tj] 


K»1 


(3a) 


f,  (*o . h)' 


— 1 - expf-trlM-r'T,)] 


K*I 


jn  detMj. 

where  the  sample  covariance  matrices  To  and  Tj  are 


(3b) 
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t  , = ^|(Iq  -  m*0  -  k/1 + 1 

The  likelihood  ratio  test  is  based  on  separate  maximizations 
of  fo  and  fj  over  the  unknown  parameters  My  and  b.  The 
detection  statistic  is  then  taken  as  the  ratio  of  the  maxima. 
The  estimates  of  My  and  b  which  maximize  £o  and  f  i  are  by 
definition  the  maximum  likelihood  estimates. 

It  is  well  known  that  the  maximum  likelihood  estimate 
of  the  covariance  matrix  is  given  by  the  sample  covariance 
matrix  [7].  Hence,  To  maximizes  fo  over  My  and  Ti 
maximizes  fi  over  My.  We  now  have 

iK*l 


Mr 


(4a) 

(4b) 


max  1,  =  (c 

M.„  1 


lK*l 


[cnf  detTj] 

T  .  (5) 

b  appears  only  in  Tj  so  for  convenience  we  form  the  K+lsl 
root  of  the  ratio  of  fj  to  fo.  termed  L(b),  and  maximize  it 
over  b  to  obtain  the  likelihood  ratio  test 

detT 
L(b)=  “ 


max  L(b)  = 


detT, 


detT„ 


min  detT 


H, 


< 

Ho 


(6a) 


(6b) 


where  Lq  is  the  threshold  parameter. 

Using  a  derivation  analogous  to  Kelly's  [1],  b  is  given 
by 

h  -i 

b  =  hS.  ‘.fo. 


s^S 


(T) 


where 


«  V  H 
S  =  >  zi 

%  1  l.  (8) 

Substitute  for  b  in  (6b)  to  obtain  the  likelihood  ratio  test 

,  V  Hl 

l  -  Z0  >  l 


Hc 
+  Z0S 


•■A 


H— -t. 


*T  (9) 

Note  that  this  test  is  of  the  same  form  as  Kelly's,  as  it 
should  be  since  the  transformation  T  is  linear.  Thus,  the 
analysis  given  by  Kelly  applies  directly  to  (9)  with  minor 
modifications  due  to  T.  A  summary  of  the  properties  of 
this  test  is  given  here;  the  reader  is  referred  to  [1]  for  details. 

The  probability  of  false  alarm  (PFA)  is  given  by 


PFA 


1 V 

ij 


-0'  (10) 
This  test  has  the  constant  false  alarm  rate  (CFAR)  property 
since  the  PFA  depends  only  on  K  and  P;  it  is  independent  of 
M  and  T.  The  probability  of  detection  (PD)  is 


PD 


L  L 

I 

k  •  1 


&KA.UUII  kSJ  U 


(11) 


where  L  =  K+l-P  and 


^(y)  =  G^ryXfr)  dr 


k-l 


(12a) 


Gk(x)  =  exp(-x)  £  — 


o  n! 

-2  L 


(12b) 

W- IE£4(,  -rfV 

L'(P-2)!  (12c) 

The  parameter  a  which  appears  in  the  expression  for  PD  is  a 
function  of  the  transformation  matrix  T 


a  =  |bj2s^Mj.1sT 
=  |bj2sHT(THMT) 


— ,H 

T  s 


(13) 

Transforming  T  by  any  nonsingular  matrix  on  the  right 
does  not  change  a.  The  un transformed  detection  problem 
considered  in  [1]  is  represented  by  setting  P  =  N. 

The  above  expressions  allow  evaluation  of  the  effect  of 
T  on  detection  performance.  Mapping  onto  a  subspace 
redoes  the  number  of  data  vectors  required  for  existence  of 
the  teo.1  The  original  test  required  K  >  N;  the  transformed 
test  requires  K  >  P.  If  PFA  is  held  constant,  then  reducing 
from  N  to  P  leads  to  a  smaller  threshold  Lq  which  tends  to 
increase  PD.  r  represents  a  Beta  distributed  (12c)  loss  factor 
which  arises  due  to  the  estimation  of  the  covariance  matrix 
[1],  The  loss  factor  reduces  the  effective  value  of  a  (see 
(12a)).  Note  that  E(r)  =  (K+2-P)/(K+l).  As  P  decreases, 
the  distribution  of  the  loss  factor  becomes  concentrated 
closer  to  the  maximum  value  of  unity,  which  tends  to 
increase  PD.  This  increase  corresponds  to  the  gain  in 
statistical  stablity  associated  with  estimation  of  the 
covariance  matrix  in  the  subspace.  Much  fewer  parameters 
(P2  vs.  N2)  must  be  estimated  given  the  same  number  of 
data  vectors. 

PD  also  increases  as  a  increases.  The  dependence  of  a 
on  T  is  illustrated  by  reexpressing  (13)  in  the  form 


i  =  jbf sHM4V^f 


.H/2 


_J_Hm  .1/2...H/2. 

T\T  M  M 


t)''thm1/2]m 


1/2 


a  =  |b|sM  |M  i\t  ivi  m  r/  i  im  jm  s  ((4) 

where  M1/2MHf2  is  the  Cholesky  factorization  of  M.  The 
term  in  brackets  is  a  projection  matrix.  Thus,  a  is 
proportional  to  the  norm  of  the  projection  of  M~,/2s  onto 
the  P  dimensional  space  spanned  by  MH/2T.  This  norm  is 
upper  bounded  by  T  nonsingular  (P  *  N).  If  M  were 
known,  we  could  attain  the  upper  bound  with  P  =  1  by 
choosing  T  =  M~ls.  However,  M  is  not  known  and  we 
will  not  in  general  attain  the  upper  bound.  This  results  in  a 
loss  in  performance  associated  with  subspace  processing. 

The  parameter  a  is  related  to  the  signal  to  noise  ratio  at 
the  output  of  a  linearly  constrained  minimum  variance 
beamformer.  Using  the  freedom  to  transform  T  on  the  right 
by  a  nonsingular  matrix,  we  assume  without  loss  of 
generality  that  T  is  of  the  form 

T  s  [$Tj  ]  (15) 

where  the  N  by  P-1  matrix  Ts  satisfies  sHT$  =  0.  The 
partitioning  in  (15)  is  always  possible  provided  sHT  *  0,  a 
condition  which  is  satisfied  by  all  T  of  interest  since  a  =  0 
if  shT  *  0.  Substituting  this  T  into  (13)  yields 


■WH* 


HSW  sHMs  sHl 
\.  T^Ms  T? 


sHMT, 
MT, 


hl\ 


□ 


1,1 


(16) 


odeju 


/or 


:ial 


r 
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where  (A)i,i  denotes  the  element  in  the  first  row  and  first 
column  of  A.  Apply  the  formula  for  the  inverse  of  a 
partitioned  matrix  to  obtain 

W2{shs)2 

d  -  1  ’  -  — 

sHMs -sHMT^MT.r^^Ms  (!7) 
a  represents  the  signal  to  noise  ratio  at  the  output  of  a 
partially  adaptive  linearly  constrained  minimum  variance 
beamformer  with  weight  vector  w  expressed  in  generalized 
sidelobe  canceller  (GSC)  [2]  form  as  w  =  s  -  Tsw„.  In 
GSC  terminology,  s  is  a  fixed  beamformer  designed  to  pass 
the  signal,  Ts  is  the  signal  blocking  matrix  which  is 
orthogonal  to  the  signal,  and  wn  is  a  set  of  unconstrained 
adaptive  weights.  wn  is  chosen  to  minimize  the  output 


power 


min (s  -  T sW,,^  {xxH  }(s  -  T swn) 


which  has  solution 

wn  =(t^MTs)_1T^Ms  (19) 

This  solution  is  obtained  independently  of  whether  x  is 

distributed  according  to  Hi  or  H0  because  of  the 
orthogonality  of  Ts  and  s.  The  output  power  due  to  the 
interference,  Pj,  is  given  by 

P,  =  sHMs  -sHMT  (tHMT  )  THMs  ™ 

1  s'  s  s'  s  (20) 

and,  assuming  the  signal  is  present  in  x,  the  output  power 
due  to  the  signal,  Ps,  is 

'.-wW.  C1) 

Thus,  (17)  is  rewriuen 


3.  Design  of  T 

The  discussion  in  the  previous  section  suggests  a 
general  philosophy  for  designing  T  to  optimize  the 
subspace  based  detector's  performance.  P  should  be  chosen 
as  small  as  possible  and  T  should  lead  to  values  of  a  which 
are  close  to  the  upper  bound.  Below  we  propose  two 
similar  methods  for  designing  T. 

Let  0  be  a  vector  which  parameterizes  the 
interference/noise  environment  and  denote  the  corresponding 
covariance  matrix  as  M(0).  0  can  represent  interferer 
locations,  power  levels,  isotropic  background  noise 
characteristics,  etc.  For  example,  if  we  assume  the 
interference  environment  consists  of  two  interferers  in  white 
noise,  then  0  would  represent  the  interferer  locations. 
Denote  a  as  a  function  of  0  as  a(0) 

a(0)  =  s,MTHM(0)T)  THs  (23) 

For  convenience  we  assume  Ibl »  1 .  The  goal  is  to  optimize 
T  over  a  set  of  likely  interference  scenarios  denoted  by  fl. 
For  example,  maximize  a  weighted  average  of  a(0)  over  Cl 


max  b(0)a(0)d0 


where  b(0)  is  a  nonnegative  weighting  function. 

The  correspondence  between  a  and  the  SNR  at  the 
output  of  a  partially  adaptive  beamformer  can  be  used  to 


pose  an  alternative  optimization  problem.  a(0)  is  inversely 
proportional  to  Pj(0)  (22)  so  maximizing  a(0)  over  T  is 
equivalent  to  minimizing  Pi(0)  over  T,  where  Tf  is 
constrained  to  be  orthogonal  to  the  signal  vector  s.  Here 
the  constraint  that  T  contain  a  component  in  the  space 
spanned  by  s  is  explicitly  enforced;  this  constraint  is 
implicit  in  maximization  of  a(0)  over  T  although  it  can  be 
made  explicit  by  constraining  T  to  be  of  the  form  given  in 
(15).  Minimizing  a  weighted  average  of  the  interference 
power  over  (2  gives  the  problem 

min  j  b(0)  P,(0)  d0 

1.  Je.o  (25) 

This  problem  is  identical  to  the  partially  adaptive 
beamformer  design  problem  posed  ir  [2).  Thus,  the 
methods  developed  in  [2,5,6]  can  be  employed  to  solve  (25). 

Note  that  (24)  and  (25)  are  not  equivalent.  At  a  single 
value  of  6,  maximizing  a(0)  over  T  is  equivalent  to 
minimizing  P](0)overTs.  However,  the  optimization 
problems  (24)  and  (25)  perform  the  maximization  and 
minimization  over  a  range  of  0.  Using  the  relationship 
between  a(0)  and  Pj(0)  we  see  that  (25)  is  equivalent  to 


b(0)-J— d0 
,  a(0) 


Thus,  while  (24)  maximizes  the  average  SNR,  (25) 
minimizes  the  average  of  the  ;nverse  SNR.  The  inverse 
SNR  criterion  more  heavily  weights  small  SNRs  and 
deemphasizes  large  SNRs  in  choosing  T. 

The  remainder  of  this  paper  considers  the  optimization 
problem  of  (24).  A  closed  form  solution  ro  (24)  is  not 
apparent  so  we  follow  a  strategy  similar  to  that  posed  for 
panially  adaptive  beamformer  design  [2],  Approximate 
solutions  are  obtained  by  designing  T  one  column  at  a  time. 
Partition  T  as  T  =  [  t  T0  ]  where  t  is  an  N  dimensional 
vector  and  T0  is  assumed  known.  The  P  columns  of  T  are 
designed  one  at  a  time,  with  each  new  column  being  a 
function  of  previously  designed  columns.  a(0)  is  expressed 
as  a  function  of  t  and  T0  using  the  relationship  for  the 
inverse  of  partitioned  matrices  as 

a(0)  =  s’VW-  cH(0)c(0)+  *  A  (H6)CH(e)c  (6)A(6)t 

t  A  (0)A(0)t  07) 


P(0)  =  I  -  MH/2(0)To(T>(0)To)  'AAo)  ^ 

A(0)  =  P(0)MH/2(0);  c(0)  *  P(0)M  1/2(0)s  (28b) 

The  first  two  terms  in  (27)  do  not  depend  on  t.  Thus,  the 
optimization  problem  is  expressed  as 


( 

nax  I 

'  L 


tV(0)c(0)cH(0)A(0)t  ^ 
tHAH(0)A(0R 


where  represents  a  subset  of  Q. 

Development  of  an  approximate  solution  to  an 
optimization  problem  of  the  same  form  as  (29)  is  given  in 
[2]  and  [5].  Application  of  these  results  leads  to  a  t  which 
satisfies  the  set  of  linear  equations 

At  =  c  (30) 

where 
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A  = 


c  = 


j  b(8)  AH(8)A(6>  d8 

J»«o, 

■ 

(31a) 

b(8)AH(8)c(8)de 

.  9*0, 

(31b) 

Consider  the  case  where  the  design  region  is 
approximated  by  a  single  point  of  the  parameter  space,  = 
90  such  that  the  integrals  in  (31a,b)  are  well  approximated 
by  the  integrands  evaluated  a:  80-  Now  A  = 

b(9o)AH(9o)A(0o)  and  c  =  b(90)AH(9o)c(9o).  The  presence 
of  the  projection  matrix  P(-i  ,)  in  the  definitions  of  A(90) 
and  c(90)  implies  that  c(90)  ies  in  the  space  spanned  by  the 
columns  of  A(90)  and  that  A  is  not  full  rank.  Let  A  = 

USVH  be  the  singular  value  decomposition  of  A.  The 
minimum  norm  solution  to  (30)  is 

*  ■  Vl"'uHC(9#)  (32) 

where  here  superscript  -1  denotes  pseudo  inverse. 

Compose  U  of  the  columns  of  U  corresponding  to 
nonzero  singular  values.  Using  (32)  in  (27)  we  have 

tHAH(9o)A(9o)t=cH(9o)UU  c(6J 


so  that 


tHA>o)c(9o)=cH(9o)UU  c(9o) 

=  iV(9ji-  c"(9j(l  -  UUH)c(9o) 


(33) 


(34) 


Now  U  is  a  basis  for  the  space  spanned  by  the  columns  of 
A(0O).  c(0o)  lies  in  this  space  so  the  projection  of  c(0o) 


onto  the  space  orthogonal  to  U  is  zero  and  (34)  simplifies 
to 

a(0o)  =  sHm  (90)s.  (35) 

which  is  the  maximum  value  for  a(0o)- 

Thus,  as  the  design  region  shrinks  to  a  point,  the 
design  procedure  suggested  here  leads  to  a  value  for  a(0o) 
which  is  equal  to  the  value  obtained  in  the  absence  of 
subspace  processing.  This  suggests  that  the  best  overall 
performance  is  obtained  by  designing  each  column  of  T  over 
separate  subregions  of  Q  which  are  are  made  as  small  as 
possible.  If  these  subregions  are  sufficiently  small,  then  the 
subspace  based  detector  should  show  very  little  degradation 
in  a(6). 

The  correspondence  between  a  and  the  SNR  at  the 
output  of  a  partially  adaptive  beamformer  facilitates  analysis 
of  the  sensitivity  of  the  design  procedure  to  differences 
between  the  actual  interference  environment  and  the  class  of 
environments  assumed  during  design.  This  allows 
determination  of  which  parameters  are  most  important  to 
include  in  8.  If  the  design  results  in  a(6)  being  large  for  6 
in  the  set  Q,  then  Pi(8)  will  be  small  for  8  in  the  set  £2. 
We  address  the  sensitivity  of  a(8)  to  environments  8  not  in 
n  by  examining  the  sensitivity  of  P](8)  to  environments  8 
is  not  in  Q. 

Suppose  the  interference  environment  consists  of  J 
point  interferes  in  white  noise  so  that 


i  *  1  (36) 


where  d(a)  is  the  array  response  or  direction  vector 
corresponding  to  direction  a.  In  general,  6  for  this 
environment  would  include  L,  <JW2,  a-;2,  and  cti,  i  =  1,2, 
....  L.  Pi(8)  is  expressed  as 

Pl(e>.(s  -  T  .JHM(8)(s  -  T  wn) 

-  .1 4Hd<“i>  -  -KW  *  4  -  T,«J2 
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Note  that  wMd(a)  is  the  response  in  direction  a  associated 
with  w. 

There  are  two  factors  in  Pj:  the  interference  term  depend 
on  how  well  the  response  of  s  matches  the  response  of 
T,wn  at  the  interferer  directions  a;,  and  the  white  noise 
term  depends  on  the  norm  of  s  -  Tswn.  In  general  these 
two  terms  are  inversely  related;  as  response  matching 
improves  the  norm  gets  larger.  The  relative  sizes  of  a,2  and 
ow2  determine  the  relative  importance  of  response  matching 
versus  weight  vector  norm. 

If  the  Cj2  are  much  larger  than  aw2,  then  response 
matching  dominates.  Given  a  Ts  which  provides  good 
response  matching.  Pi  will  be  insensitive  to  variations  in 
o,2.  Pi  will  be  sensitive  to  variations  in  the  interferer 
directions.  If  Ts  is  chosen  to  provide  good  response 
matching  over  one  range  of  a,  it  is  unlikely  that  good 
response  matching  will  be  obtained  over  a  different  range  of 
a.  Pi  should  not  be  sensitive  to  changes  in  the  number  of 
interferers  L  provided  L  <  P-1.  Assuming  the  d(a,),  i  * 
1.2.....L  are  linearly  independent,  the  response  matching 
term  represents  L  equations  in  P-1  unknowns  (»„).  If  L  < 
P-I,  then  the  response  matching  term  can  always  be  made 
equal  to  zero  independent  of  Ts,  although  this  may  result  in 
a  large  norm  for  s  -T,wn.  A  good  Ts  will  achieve 
response  matching  close  to  zero  while  keeping  the  norm 
small. 

If  ow2  is  much  larger  than  the  a;2,  then  the  norm 
dominates  and  the  design  will  not  be  sensitive  to  variations 
in  a;2,  to  L,  or  the  interferer  locations.  A  Ts  designed  for 
Oj2  much  larger  than  aw2  should  perform  well  when  aw2  is 
much  larger  than  the  O;2  because  a  Ts  which  results  in  good 
response  matching  can  yield  a  small  norm  by  setting  wn  = 
0.  In  contrast,  a  Ts  designed  for  ow2  much  larger  than  the 
Oj2  is  not  necessarily  capable  of  good  response  matching. 

The  preceding  discussion  suggests  that  the  most 
important  parameters  to  include  in  6  are  the  interferer 
directions.  The  best  all  around  performance  is  obtained  by 
designing  T,  assuming  the  Oj2  is  much  larger  than  ow2. 
The  exact  number  of  interferers  does  not  appear  to  be 
significant  as  long  as  L  <  P-1. 

4.  Simulations 

A  linear  equal  spaced  array  of  30  sensors  is  used  to 
illustrate  the  potential  performance  improvements  which  are 
possible  by  paforming  detection  in  subspaces.  The  sensors 
are  spaced  at  one  half  wavelength.  The  signal  to  be  detected 
arrives  from  the  direction  perpendicular  to  the  array  so 
s  *  SO-1/2!  1  1  1  ...  1]H. 
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The  magnitude  of  b  is  set  to  unity.  The  noise  consists  of 
two  interferers  in  white  noise,  where  the  ratio  of  interferer 
power  to  white  noise  power  is  30  dB.  This  implies 

M  =  1000d(a1)dH(a1)  +  lOOGJfa^fCj)  +  50_l/2I  (3g) 

where  d(a)  is  assumed  to  have  unit  norm  and  aj,  02 
represent  the  interferer  directions 

The  subspace  based  detector  is  designed  for  P  =  4  by 
solving  (30)  for  each  of  the  columns  of  T.  b(0)  is  equal  to 
unity  and  the  integrals  in  (3la,b)  are  approximated  by  sums. 
The  elements  of  6  are  a  1  and  02  and  T  is  designed  to 
optimize  performance  over  the  range  .1  <  sinaj  <  1,  .1  <, 
sina2  <  1.  The  subregion  over  which  each  column  of  T  is 
designed  is  given  in  Table  1.  These  regions  were  not 
customized  to  optimize  performance,  but  were  selected  in  a 
somewhat  arbitrary  manner. 

The  values  of  a  for  the  subspace  and  non-subspace 
detectors  are  computed  at  182  different  combinations  of  ai 
and  a;  or.  the  region  .1  <  sinai  <  1,  .1  <  sina2  £  1. 
Table  2  presents  a  histogram  of  the  ratio  of  non- subspace  a 
to  subspace  a  measured  in  dB.  This  number  is  a  measure  of 
the  loss  in  a  resulting  from  subspace  processing.  The 
design  procedure  is  clearly  effective  since  177  of  the  182 
cases  show  less  than  1  dB  loss  due  to  the  mapping  from  a 
50  dimensional  space  to  a  4  dimensional  space.  The  worst 
case  loss  is  3.51  dB. 

The  PD's  of  both  detectors  are  computed  for  the  same 
182  cases  assuming  K  =  75  signal  free  data  vectors  and  a 
constant  PFA  =  10-6.  The  PD's  for  the  non-subspace 
detector  ranged  from  a  minimum  of  .4855  to  a  maximum  of 
.5066.  Table  2  presents  a  histogram  of  the  PD's  for  the 
subspace  detector.  Subspace  processing  results  in  a 
tremendous  improvement  in  the  PD;  the  PD  is  greater  than 
.999  for  177  of  the  182  cases.  The  worst  case  PD  of  .8539 
is  significantly  greater  than  the  best  case  non-subspace  PD 
of  .5066.  The  two  cases  where  the  PD  is  less  than  .9 
correspond  to  the  two  cases  where  the  loss  in  a  is  greater 
than  1.5  dB.  Customization  of  the  regions  used  to  design 
each  column  of  T  could  be  used  to  improve  a  and  the  PD  at 
these  points. 

A  histogram  of  the  PD's  over  the  same  set  of 
interference  scenarios  is  given  in  Table  4  assuming  PFA  = 
10-6  and  K  =  30.  The  non-subspace  detector  does  not  exist 
in  this  case  because  K  <  N.  The  worst  case  PD  is  .6472; 
175  of  the  cases  have  PD's  better  than  .99.  The  subspace 
detector  performs  significantly  better  than  the  non-subspace 
detector  even  though  it  has  less  than  one  half  the  number  of 
signal  free  data  vectors  available. 

5.  Summary 

The  detection  problem  posed  by  Kelly  [1  ]  is  modified  to 
perform  detection  on  data  which  is  mapped  into  a  subspace 
prior  to  processing.  The  detection  performance  tends  to 
increase  due  to  the  reduction  in  data  dimension,  but  tends  to 
decrease  due  to  a  loss  in  SNR  associated  with  the  mapping 
into  the  subspace.  A  procedure  is  proposed  for  designing 
the  subspace  transformation  to  minimize  the  SNR  loss. 
The  subspace  design  problem  for  optimizing  detection  is 
shown  to  be  closely  related  to  the  partially  adaptive 
beam  former  design  problem.  Simulations  illustrate  the 
effectiveness  of  subspace  detection.  The  gain  in  detection 
performance  associated  with  reducing  dimension  far  exceeds 


the  detection  loss  associated  with  the  loss  of  SNR  in  the 
subspace  for  tbe  cases  studied  here. 
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Table  1.  Design  regions  for  the  columns  of  T. 


Column 

sinai 

sina2 

1 

.6  -  1.0 

.6  -  1.0 

2 

.6  -  1.0 

.1  -  .58 

3 

.3  -  .58 

.1  -  .58 

4 

.1  -  .28 

.1  -  .28 

Table  2.  Difference  between  101og(a)  for  the  non- 
subspace  detector  and  loiog(a)  for  the  subspace  detector. 


Range  for  Difference 

Number  of  Cases 

0  <  dif  <  .5 

174 

.5  <  dif  <  1.0 

3 

1.0  <  dif  <  1.5 

3 

1.5  <  dif  <3.6 

2 

Table  3.  PD's  for  subspace  detector  with  K  =  75.  The 
non-subspace  detector  PD's  range  from  .4855  to  .5066- 


Range  for  PD 

Number  of  Cases 

.85  <  PD  <  .9 

2 

.9  <  PD  <  .99 

0 

.99  <  PD  <  .999 

3 

.999  <  PD  <  .9999 

4 

.9999  <  PD  <  .99999 

173 

Table  4.  PD's  for  subspace  detector  with  K  =  30. 


Range  for  PD 

Number  of  Cases 

.64  <  PD  <  .9 

2 

.9  <  PD  <  .99 

4 

.99  <  PD  <  .999 

175 

167 


